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ABSTRACT 



In this paper we propose a new plausible mechanism of supernova explosions spe- 
cific to close binary systems. The starting point is the common envelope phase in 
the evolution of a binary consisting of a red super giant and a neutron star. As the 
neutron star spirals towards the center of its companion it spins up via disk accretion. 
Depending on the specific angular momentum of gas captured by the neutron star via 
the Bondi-Hoyle mechanism, it may reach millisecond periods either when it is still 
inside the common envelope or after it has merged with the companion core. The high 
accretion rate may result in strong differential rotation of the neutron star and gener- 
ation of the magnetar-strength magnetic field. The magnetar wind can blow away the 
common envelope if its magnetic field is as strong as 10^^ G, and can destroy the entire 
companion if it is as strong as 10^^ G. The total explosion energy can be comparable 
to the rotational energy of a millisecond pulsar and reach 10^^ erg. However, only a 
small amount of ^^Ni is expected to be produced this way. The result is an unusual 
type-II supernova with very high luminosity during the plateau phase, followed by a 
sharp drop in brightness and a steep light-curve tail. The remnant is either a solitary 
magnetar or a close binary involving a Wolf-Rayet star and a magnetar. When this 
Wolf-Rayet star explodes this will be a third supernovae explosion in the same binary. 



INTRODUCTION 



Usually, neutron stars (NS) have magnetic field B ^ 10^^ G 
and rotate with a period of fraction of a second, the Crab 
pulsar being a typical example. However, we now know that 
the "zoo" of NS is much more diverse and they can have 
both much weaker and much stronger magnetic field and 
rotate with both much longer and much shorter periods. 

Around 10% of NS have surface di polar magnetic field 
B - lO^'^ - lO^^G dKouveliotou et al.|[l998J . These "mag- 
netars" are believed to be born in core collapse explo- 
sions of rapidly rotating stars. During the collapse, the 
proto-neutron star naturally develops strong differential ro- 
tation, which allows for generation of magnetic field via a-Q- 
dynamo ( ,Duncan Thompson|1992 Thompson & Duncan 



e.g. Usov| 199*2 Thompson et al.|2004| . Turbulence required 
for the magnetic dynamo action can also be generated via 



the magneto-rotational instability (MRI, Balbus & Hawley 



1991 ). Calculations based on the linear theory show that in 



the supernova context strong saturation field can be reached 
very quickly, on the time scale of only several tens of rota- 
tional periods (Akiyama et al. 2003| Obergaulinger et al. 
2009). 

Millisecond pulsars are found in low mass binaries, and 
it is generally thought that they have been spun up via disk 



accretion (Alpar et al. 1982" Archibald et al. 2009). This 



1993). The strength of saturated magnetic field strongly de- 



pends on the rotational period, with shorter periods leading 
to stronger magnetic field and more rapid release of rota- 
tional energy. In order to generate the magnetar strength 
magnetic field the rotation period must be around few mil- 
liseconds ( [Duncan Thompson p!992| ). Such a strong field 
allows to release up to 10 erg of magnetar 's rotational en- 
ergy in very short period of time. This is sufiicient to drive 
extremely powerful supernova explosions, on the hypernova 
scale, and to produce powerful Gamma Ray Bursts (GRB, 



origin implies mass increase by about O.2M0 compared to 
normal radio pulsars, whose masses are narrowly distributed 
around 1.35Mq (T horsett Chakrabarty|1999t . It is rather 
difiicult to measure the mass of a millisecond pulsar, but the 
few available results agree with this prediction of the accre- 
tion model (Kaspi et al.|1994 Jacoby et al.|2005 Demorest 
et al.|2010| ). The most massive pulsar found to date, almost 



2M0, is a millisecond pulsar (Demorest et al. 2010). The 



magnetic field of these millisecond pulsars is very low, down 
to 10^ Gauss. Most likely, their initial magnetic field was 
of similar strength to normal pulsars, but now it is buried 
under the layers of accreted matter ( [Bisnovatyi-Kogan 
Romberg 1974 Alpar et al. 1982). The reason why these 
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pulsars could not generate magnetar-strength magnetic field 
in the same way as in the core-collapse scenario is the very 
long time scale of spinning up compared to the viscous time- 
scale. As the result, the differential rotation remains weak 
and there is not enough energy for effective magnetic dy- 
namo ( |Spruit|1999 ). 

The rotational frequency of the fastest X-ray pulsar is 
~ 760 Hz (iChakrabarty et al. 2003), whereas the fastest 
known radio pulsar has the frequency of 641 Hz (Backer 



et al.||T982 ). Most likely, it is the gravitational radiation 
losses what places the upper limit on the rotation rates be- 
cause at high spin the r-mode oscillations become excited 
( Shapiro T eukolsky| [l983| |Levin||1999t . [Spruit (1999) ar- 
gued that this instability may also result in magnetic ex- 
plosion. The idea is that the heating of NS, associated with 
these oscillations, reduces its viscosity and decouples its inte- 
rior from the outer layers. Being most disturbed, the outer 
layers rapidly loose some of their angular momentum via 
gravitational radiation. This leads to strong differential ro- 
tation and generation of magnetar strength magnetic field 
in the NS interior. This field becomes unstable to buoy- 
ancy, emerges on the surface, and magnetically driven pul- 
sar wind rapidly extracts the rotational energy of the NS. 
Spruit] (|^1999) proposed this as an alternative scenario for 
long GRBs. It is unlikely that a supernova- like event can 
accompany a GRB in this scenario. Although the wind en- 
ergetics is sufficient, only a small fraction of this energy can 
be deposited into the companion star, simply because of 
its small geometrical cross-section. Moreover, recent results 
suggest that the amplitude of r-modes may saturate at a 
much lower level due to nonlinear interaction with other 
modes (Arras et al.|2003| [Br ink et al. 2005; Bondarescu et al 
2007|. — — 



A similar recycling of NS may occur during the common 
envelop (CE) phase, after the primary becomes a red super 
giant (R SG; [Bisnova tyi-Koga n Syunyaev[[1971[ [Paczyn- 
ski| 19761 [Tutukov Yungelson|1979| [Postnov Yimgelson 
2006). Due to the dynamic friction, the NS then spirals in- 
side the RSG, accreting on its way. Now one can imagine 
two interesting outcomes of such process. First, the neu- 
tron star may accumulate too much mass and collapse into 
a black hole (BH). This BH is likely to be rapidly rotating 



and drive a stellar explosion in the collapsar fashion (Fryer 



&: Woosley|1998l [Zhang Fryer|200r[ [Barkov Komissarov 

2oTof ^ 

Second, the NS may first spin up to a millisecond pe- 
riod and drive a magnetic explosion of the type proposed 



by Spruit (1999) but now inside the common envelope. The 



magnetar wind can keep energising such supernovae, pro- 
ducing a similar effect to radioactive decay ( Woosley] [2010 



Kasen Bildsten|[2010 ). High accretion rates may modify 
the way the NS is recycled. The accreted gas can form a mas- 



sive rapidly rotating layer above the NS crust (Inogamov & 
Sunyaev 1999[ 2010). The strong differential rotation be- 
tween the layer and the NS core may result in development 
of the Kelvin- Helmholtz instability when the NS crust melts 
down under the weight of the layer. This may lead to tur- 
bulence and strong amplification of the NS magnetic field. 
The accretion onto neutron stars during the in-spiral 



keeps the mass accreting rate well below the rate of the 
Bondi-Hoyle capture. As the result, the NS accumulates very 
little mass while still inside the common envelope. On the 
other hand, he suggested that during the merger with the 
companion core the mass accretion rate rises and the NS 
collapses into a black hole. While carefully analysing var- 



has been studied by Chevalier ( 1996 J , who concluded that 
the high angular momentum of the gas gravitationally cap- 
tured by NS prevents it from effective neutrino cooling and 



ious effects of rotation. Chevalier (1996) did not use any 
particular model for the primary. Moreover, there is a great 
deal of uncertainty with the regard to the specific angular 
momentum of the gravitationally captured gas. In our study 
we come back to this problem, consider realistic models of 
RSG stars and allow for the uncertainty. 



The paper is organised as follows. In Section 2 we con- 
sider the accretion and recycling of NSs during the in-spiral. 
We conclude that the outcome is very sensitive to the as- 
sumed specific angular momentum of the gravitationally 
captured gas. Given the current uncertainty with the regard 
to the angular momentum, it seems possible that NSs be- 
gin to accrete with Bondi-Hoyle rate while still inside com- 
mon envelopes. In this case, they rapidly spin up to mil- 
lisecond periods. In Section 3 we speculate on the possible 
mechanisms of generating magnetar-strength magnetic field 
and study the magnetic interaction of millisecond magnetars 
with accretion flows typical for the in-spiral problem. We 
find that if the magnetar forms inside the common envelope 
and its magnetic field is about 10^^ G, the magnetospheric 
pressure can overcome the gravity and drive an outflow, with 
eventual release of up to lO^^erg of the magnetar rotational 
energy. If the magnetar forms only after the merger with 
the core then higher magnetic field, ^ 10^^ G, is required to 
drive stellar explosions. In Section 4 we discuss the prop- 
erties of such explosions and their observational signatures. 
Because the RSG envelope is rich in hydrogen the supernova 
will be classified as type-II. Because of the very high en- 
ergy release and relatively small mass of the RSG, the speed 
of the ejecta is expected to be very high, ^ lO^cm/s, and 
the luminosity at the plateau phase ~ lO^^erg/s. However, 
due to the small amount of ^^Ni generated in the explosion 
and the rapid spin-down of the magnetar, the plateau is fol- 
lowed by a sharp drop in brightness and a steep light-curve 
tail. We show that the termination shock of the magnetar 
wind produces a high energy synchrotron and inverse Comp- 
ton emission which may energise the tail, but the supernova 
ejecta soon becomes transparent to this emission, limiting 
its potential to mimic the effect of radioactive decay. The 
flux of gamma-ray emission is expected to be rather low 
and difficult to observe, unless the explosion occurs in the 
Local Group of galaxies. In the case of off-center explosion, 
the remnant is a very close binary system consisting of a 
WR star and a magnetar, but the strong magnetic field pre- 
vents magnetar from accreting plasma of the WR wind. Our 
results are summarized in Section 5. 



In this paper, the dimensional estimates are presented 
using the following notation: the time tx,n is measured in 
lO'^s, the distance i?x,n in 10"^ cm, the speed K,n in 10"^ cm/s, 
the mass Mx,n in 10^ M©, mass accretion rate Mx,n in 
WMq/jt, and the magnetic field Bx,n in lO^'G. 
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2 IN-SPIRAL DYNAMICS AND RECYCLING 
OF NEUTRON STARS 

As the neutron spirals inside its giant companion it accretes 
mass and angular momentum. The accretion can proceed 
at the Eddington rate or at the much higher Bondi-Hoyle 
rate. This depends on whether the radiation is trapped in 
the accretion flow and the neutrino cooling is sufficiently 



effective (Houck & Chevalier 1991 Chevalier 1993 1996 



Fryer et al. 1996). Only if the accretion proceeds at the 



Bondi-Hoyle rate the recycling of NS is sufficiently fast and 
the explosion can occur during the inspiral. 

In the Bondi-Hoyle problem, an accreting point mass M 
is moving with finite speed Voo through a uniform medium 
of mass density poo and sound speed Cs,oo- The Bondi-Hoyle 
mass accretion rate is well approximated by the equation 



Mbh - TtRaPooVo 



2GM 



where 

is the accretion radius, and 

^2 \ 3/4 



where A4c 



Mi 



1 + 



(1) 



(2) 



(3) 



d/cs,oo is the Mach number (Shima et al. 



1985). The mean angular momentum of accreted matter is 



The in-spiral problem is more complicated due to the 
density and velocity gradients across the direction of motion 
of the NS. Provided, the gradients are small on the scale of 
accretion radius the above expression for M is still quite ac- 
curate pshii et al.|[r993 i Ruffert 1997 , 1999 ). However, the 
accreted matter can now have non- vanishing mean angular 
momentum. This was first suggested by |Illarionov SuiT] 
yaev (1975) in connection with X-ray binaries, where the 
accretion disk of the black hole is fed by the stellar wind 
of the primary. Taking into account only the velocity gra- 
dients related to the orbital motion, they found that the 
mean specific angular momentum of accreted matter inside 
the accretion cylinder is 



(4) 



where Q, is the angular velocity of the orbital motion. 



Shapiro & Lightman (1976) took into account the density 



gradient of the wind as well and refined this result. 



(5) 



where the sign minus shows that the disk rotation is now 
in the opposite sense to the orbital motion. Later, [Pavies fc| 
Pringle ( 1980 ) criticised the approach used by the previous 



authors. They generalised the original calculations of |Hoyle] 
Lyttletoni (^1939) to include the density and velocity gra- 
dients and concluded that to the first order in RA/h, where 
h is the length scale of the variations, (Ja) = 0. The 2D 
numerical simulations of |Ishii et al.| (p!993| seemed to agree 
with this conclusion, whereas their 3D simulations still in- 



dicated net accretion of angular momentum. Ruffert j ( 1997 



1999) carried out extensive 3D numerical investigation of 



the problem. Introducing angle-dependent accretion radius. 



he also derived another version of the analytic expression for 



(Ja) = ^(6ev + ep)voRA, 
where 



(6) 



(7) 



Ra dp Ra dv 

^ p dx^ ^ V dx 

are the dimensionless gradients of density and velocity re- 
spectively, with X being the transverse Cartesian coordinate. 
However, the numerical results did not agree with this ex- 
pression, with (Ja) varying between 7% and 70% of what is 
available in the accretion cylinder. Moreover, in the simula- 
tions only the cases with either velocity or density gradient 
were considered, whereas in the in- spiral problem both are 
present. Although the computational resources and numer- 
ical algorithms have improved since 1999, no further attack 
on this problem has been attempted yet and the issue of 
(Ja) remains open. 

For the inspiral problem vq = Qa, v(x) = vq -\- Qx, and 
Cv = Ra/cl, where a is the orbital radius of NS. Assuming 
power law for density distribution inside the RSG, p oc r~^, 
we also find Cp = —jiRa/cl- This allows us to write Eq|6]as 



Oa) = -(6-n)0i?i, 



(8) 



which has the same dependence on Q and Ra as in Eqs[4]and 
|5] Given the unsettled nature of this issue, we will assume 
that 



Oa) 



(9) 



where 77 is a free parameter, which refiects our current igno- 
rance. 

In our case, M is the mass of the NS and Voo is its orbital 
speed. Assuming that the interaction does not disturb the 
RSG inside the orbit, we have 



GM./a, 



(10) 



where a is the orbital radius and is the mass of RSG 
inside the orbit. Under the same assumption, poo and Cs,oo 
can be replaced with the density and sound speed of the 
undisturbed RSG at radius r — a. The motion of the NS 
inside RSG is only mildly supersonic, ~ 1.4 (Chevalier 



1993), allowing us to write 



(11) 



and only weakly depends on the model of 



where /3 :^ O.l 
RSG. 

Given the specific angular momentum we can estimate 
the distance from the NS at which the Bondi-Hoyle trapped 
gas will form an accretion disk as 



Rc 



Ml 



2 /o4 

: arj p 



Mns 



(12) 



In the literature this radius is often called the circularisation 
radius. Since the NS is not a point mass the accretion disk 
forms only if Rc > Rns- When this condition is satisfied, the 
recycling of NS proceeds at the highest rate. 

For R ^ Rc the accretion on NS may proceed in more 
or less spherical fashion and is well described by the Bondi 
solution, in which the accretion fiow becomes supersonic at 
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Race = 0.25i?A (for the polytropic index 7 = 4/3). Its colli- 
sion with the NS and/or its accretion disk creates a shock 
wave which re-heats the flow. This will have no effect on the 
mass accretion rate only if the following two conditions are 
satisfied (Chevalier"" 1996"). 

Firstly, the radiation produced by the gas heated at the 
accretion shock should not be able to escape beyond Ra- If 
it does the mass accretion rate is limited by the Eddington 
value 



Me 



' 3xlO"^M0yr" 



This condition is satisfied when the accretion shock radius 
i?sh ~ 3 X 10^ RlfM-^-^^ cm (13) 
is smaller compared to the radiation trapping radius 

(14) 



Rtr ~ 3.4xlO^^Mocm 



For this to occur Mbh has to exceed the critical value 
Mcni ^ 1.1 X 10-^ rI^Mq yr"^ (15) 



(Chevalier 1996| ). This result was derived assuming that 
Rc > Rns and ignoring the General Relativity corrections. 

Secondly, the shock has to be inside the sonic surface of 
the Bondi flo^^ Otherwise, this shock cannot be a part of 
the stationary solution and hence the Bondi-type accretion 
cannot be realised. This condition is satisfied when the Mbh 
exceeds the critical value 



M. 



io^r: 



2.7 o4 



,,Mq yr-^ (16) 

( Chevalier|p96l ). As we shell see this condition is more re- 
strictive. 

A note of caution has to be made at this point. The 



analytical expression for the accretion shock radius ( 13 ) is 
based on a number of assumptions and simplifications and 
its accuracy has not yet been tested via numerical simula- 
tions. At the moment, this can only be considered as an 
order of magnitude estimate. When R^ becomes as small as 
the NS radius, this expression gives i?sh which exceeds the 
value given by the more reliable spherically symmetric model 



of Houck Chevalier] ( |199 1 ) by a factor of ten (see Eqs. 
and 23). Two reasons explain this disagreement. First 



13 
in 

the spherically symmetric models the gravity is relativistic, 
whereas in the models with rotation it is Newtonian. Sec- 
ond, the models use different assumptions on the geometry 
of the neutrino cooling region. It appears that Eq[T3] be- 
comes increasingly less accurate when Rc approaches i?Ns, 
overestimating i?sh- As the result, Eq |16| overestimates the 
value of Mcn2- We also note, that the critical rates are quite 
sensitive to the parameter 77, 



Mcnl oc ?7 



and M, 



cn2 OC ?7 



In this study, we have considered six different models 
of RSG. In the first two models we assumed that mass dis- 
tribution in the stellar envelope is the power law , p oc r~^, 
with n = 2 (model RSG2) and n = 3 (model RSG3). In both 
these models, the total stellar mass M* = 25M0, the helium 
core mass Mc = 1 / 3M* , the inner and the outer radii of the 



power law envelope are 7x10^ cm and 7 x 10^^ cm respec- 
tively. For these models all calculations can be made analyt- 
ically. Not only this helps to develop a feel for the problem, 
but also provides useful test cases for numerical subroutines. 
The next two models are based on stellar evolution calcu- 
lations and describe RSG in the middle of the He-burning 
phase ( |Heger et ar]|2QQ4| [Woosley et al.||2002| . Their Zero 
Age Main Sequence (AMS) masses are M = 2OM0 (model 
RSG20He) and M = 25M0 (model RSG25He). These data 
were kindly provided to us by Alexander Heger. The last two 
models, RSG20 and RSG25, describe pre-supernova RSG 
with the same AMS masse^ Their have very extended en- 
velopes, increasing the chance of CE phase. 

Figure [1] shows the Bondi-Hoyle and the critical accre- 
tion rates for these six models and the two extreme values of 
the parameter 77, representing the cases of both very efficient 
(?7 = 1) and rather inefficient (77 = 0.1) accretion of angu- 
lar momentum (see Eq|9]). One can see that the condition 
of radiation trapping is not very restrictive. For 77 = 1 it is 
satisfied for R < 3x 10^^ cm in models RSG20 and RSG25, 
for R < 3x10^^ cm in models RSG20He and RSG25He, and 
much earlier for 77 = 0.1. 

On the contrary, the condition Mbh > Mcn2 can be 
quite restrictive. In fact, it is never satisfied when 77 = 1. In 
this case, the circularisation radius is always too far away 
from the NS and the neutrino cooling is not efficient enough. 
For 77 = 0.1 the results seem more promising as this condi- 
tion is satisfied long way before NS merges with the core. 
For RSGs in the middle of the He-burning phase this oc- 
curs almost the the same point where the radiation becomes 
trapped, and a bit later for the presupernova RSGs. Thus, 
for 77 = 0.1 NS begins to hyper-accrete already inside the en- 
velope of the RSG primary. However, for such a small value 
of 77 the accretion disk becomes very compact and may even 
disappear, in which case the recycling efficiency is reduced. 

Figure|2]shows all the characteristic radii of the in-spiral 
problem for 77 = 0.1. One can see that for RSG20 and RSG25 
models, the disk accretion domain has a > 3 x 10^*^ cm. For 
the model RSG20He this domain includes the whole stellar 
envelope and for RSG25He it splits into two zones. What 
is most important is that in all models the accretion begins 
to proceed with the Bondi-Hoyle rate while still in the disk 
regime. 

Once the accretion disk is formed the specific angular 
momentum of the gas settling on the NS surface is that of 
the last Keplerian orbit. In this case the star will reach the 
rotational rate Q after accumulating 



AM : 



;0.18Moi?Ns,6^-3M( 



^0 



(17) 



where / ^ (2/5)Mi?Ns is the final moment of inertia of NS 
and Jk is the Keplerian angular momentum at the NS radius. 
Obviously, the same mass is required to recycle the observed 
millisecond pulsars in low mass binaries and this allows us 
to conclude that the NS may well reach the gravitationally 
unstable rotation rate before collapsing into a black hole. 

When the accretion proceeds at the Bondi-Hoyle rate, 
the NS mass and orbit are related approximately as 



The shock radius in Eq |l3| is obtained assuming cold supersonic ^ These models can be downloaded from the website 

flow at infinity [http://homepages.spa.umn.edu/ ^alex/stellarevolution/ ^ . 
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Figure 1. Accretion rate as a function of position inside the RSG primary for six different models: RSG2 (left top panel), RSG3 (right 
top panel), RSG20He (left middle panel), RSG25He (right middle panel), RSG20 (left bottom panel), and RSG25 (right bottom panel). 
The Bondi-Hoyle rate, Mb-h, is shown by thick solid lines, the critical value Mcni by dotted (77 = 0.1) and dashed (77 = 1) lines, and the 
critical value Mcn2 by thin solid {rj = 0.1) and dot-dashed (77 = 1) lines. 



ttinit V Mend 



(18) 



Chevalier 



( 1993 ) found that for lOM^ 



where a = —Ma/ Ma. 
RSG (7 varies slowly between 5 and 7. According to these re- 
sults, NS increases its mass by O.2M0 after its orbital radius 
decreases by only a factor of few. This shows that the "win- 



dow of opportunity" , where the disk accretion can proceed 
at the Bondi-Hoyle rate, does not have to be particularly 
wide. In order to obtain more accurate estimate, one could 
integrate the evolution equation 



AnG'^Mp 



M 



M*(l + yW-2)3/2 



,(19) 
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Figure 2. The characteristic radii of the in-spiral problem as functions of the position inside the RSG primary for the same models as in 
Figure^ and ij = 0.1. The orbital radius is shown by the solid line, the sonic point, -Race, by the dotted line, the accretion shock radius 
by the dashed line, and the circularisation radius, Rc by the dash-doted line. 



where p is the average density of the primary inside r = a, 
Cd ~ 6.5 is the drag coefficieni Shima et al. ( |1985 ), and 
^ = 1 + Mat /M* f Chevali er 1 1993t . Figures |3|and]4piow the 
results of such integration for the models RSG20, RSG25, 
and RSG25He, assuming that the hyper-accretion regime 
begins near the stellar surface and at a = 3x10^° cm respec- 
tively. One can see that in all these cases the gravitationally 



unstable rotation rate is reached when a decreases by less 
than a factor of five, in good agreement with the original 
estimate. 

For ?7 = 1, the mass accretion rate is low and the NS 
cannot spin-up significantly inside the CE. What occurs in 
this case depends on the details of CE evolution. The CE 
can either be ejected leaving behind a close NS-WR bi- 
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nary or survive, leading to a merger of the NS and the 
RSG core ( Taam Sandquist|2000| . The configuration after 
the merger is similar to that of the Thorne-Zydkow object 
(Thorne & Zytkow||1977l iBarkov et aL||2001|). The common 



view is that such objects can not exist because the effec- 
tive neutrino cooling leads to hyper-accretion and prompt 
collapse of the NS into a BH. This conclusion is based on 
the spherically symmetric model, where the circularisation 
radius is zero and the flow stagnation surface coincides with 
the surface of the NS. However, during the inspiral the or- 
bital angular momentum of the NS is converted into the spin 
of the primary. Thus, one would expect the core to be tidally 
disrupted and form a massive accretion disk around the NS. 
As the result, the cooling rate will be lower compared to 
that found in the spherically symmetric case. It is not ob- 
vious for how long this configuration may exist. However, 
this again opens the possibility of recycling of the NS up to 
the gravitationally unstable rotation rate and magnetically 
driven stellar explosion. 

In the case of core-collapse, even the magnetic field 
of magnetar strength cannot prevent accretion and fails to 
drive stellar explosion (Komissarov & Barkov 2007). The 



magnetosphere becomes locked under the pile of accret- 
ing matter, and a separated mechanism, e.g. the standard 
neutrino-driven explosion, is required to release it. However, 
the density around the Fe core of presupernovae is very high, 
leading to very high accretion rates. In our case the accre- 
tion rates are expected to be significantly lower. In the next 
section we explore if they are low enough to allow purely 
magnetically-driven explosions. 



3 MAGNETIC EXPLOSION 
3.1 Magnetic field generation 

The main route to formation of "magnetars" is believed to 
be the core collapse of rapidly rotating stars. During the col- 
lapse, the proto-neutron star naturally develops strong dif- 
ferential rotation, with the angular velocity decreasing out- 
wards. This allows generation of super-strong magnetic field 
via the a-O-dynamo ( Duncan Thompson] 1992[ Thompson 
& Duncan 1993) or magneto-rotational instability (Burrows 
et al.||2007 ). 



In our case, the conditions are rather different. The neu- 
tron star is already fully-developed, with a solid crust, when 
it enters the envelope of its companion. The typical NS crust 
has the thickness /i ~ 10^ c m and the density a t the bottom 



on temperature. When the density reaches t he critical valu e 
per ^ 10^^ g cm~^, the crust begins to melt 



Brown 



(2000), 



and the barrier separating the NS core and its outer rotating 
layer disappears. If we ignore the centrifugal force than this 
occurs when the NS accumulates mass comparable to the 
crust mass, ^ O.OSM©. The centrifugal force significantly 
reduces the effective gravitation acceleration and the above 
estimate for the layer mass is only a lower limit. Assuming 
Vcj) ~ vk ^ 10^°cm s~^ the rotational energy of the layer at 
the time of melting may well reach ~ 10^^ ergs. Once the 
crust has melted, the core and the layer can begin to inter- 
act hydromagnetically, with up to ^ ^ 10^^ ergs of energy 
in differential rotation to be utilised. 

One possibility is viscous heating and neutrino emis- 
sion. We can estimate the neutrino luminosity as L^y ^ E/td, 
where td is the neutrino diffusion time. Following the anal- 
ysis of |Lattimer &: Prakash] ( |2007t and [Barkov] ( |2008| ), this 
time can be estimated as 



' hru/c, 



(20) 



where = hpHu/Trip is the optical depth of the hot envelope 
for the neutrino emission. Assuming equilibrium with the 
electron fraction Ye — 0.05, the mean neutrino cross-section 
is about T^u ~ aueukTe/mlc^ , where ai, = 1.7 x 10~^^ cm^ 



Bahcall (1964) and ^ 3.15kTe is mean energy of neutrino 
^ Thompson et al.|2001 ). The temperature can be estimated 
as 



Te ^ (E/47varhRls) 



1/4 



-^52 '''5 -^NS,6 



where ar = 7.564 x 10 erg cm ^ is the radiation con- 
stant. The corresponding diffusion time is quite short, 

td ^ 0.1 /l5^^Pl4^52^/i?NS,6 S. 

Another process is development of the Kelvin- 
Helmholtz instability, production of turbulence and ampli- 
fication of magnetic field. For the Kelvin-Helmholtz insta- 
bility to develop the Richardson number J has to be below 
than 1/4 ( Chandrasekhar||1961 ). In our case 



J; 



: 0.2i?NS,6/^5 



(21) 



where g ~ GMns /Rns is the gravitational acceleration and 
we use Mns = 1.5M0. Thus, the instability condition is 
marginally satisfied. Once the turbulence has developed the 
magnetic field amplification is expected to proceed as dis- 



cussed in Duncan & Thompson (1992). The e- folding time 



p ^ 10 g cm (see e.g. Gnedin et al 



2001 



Lattimer & 



is given by the eddy turn-over time, which one would ex- 
pect to be below the rotation period of the outer layer. 

. ~ , 1 ■■ 1 1 1 Thus, strong magnetic field can be generated on the time 

Prakash||2007| ) . The neutron star may have already been ^^j^^ ^j^^ ^-^^^^^ neutrino diffusion time scal(|!]This 

conclusion is supported by the recent numerical simulations 
of similar problems, which show that the time as short as 
0.03 s can be sufficient to generate dynamically strong mag- 
netic field ( Obergaulinger et al.||2009| Rezzolla et al.||2011 ). 
The strong toroidal magnetic fie ld can suppers the ins tabil- 
ity in the case if pv'^ < B'^/Att d Chandr asekliarll 1 96 1 1 ) or if 



spinned up via accretion from the stellar wind, but but its 
rotation rate may still be well below the critical for develop- 
ment of the gravitational wave instability. In any case, the 
accreted gas of companion's envelope forms a very rapidly 
rotating layer above the crust. As long as the crust exists, 
this layer and the NS core are basically decoupled, and inter- 
act mainly gravitationally. A strong tangential discontinuity 
exists at the outer boundary of the crust ( [Inogamov Sun^ 
yaev|[T9 ^9 , 201Cr). 

As the mass of the outer layer increases, the crust pres- 
sure goes up. Because of the degeneracy of the crust matter, 
its pressure depends mainly on its density, and only weakly 



S > 3 X 10 G. This is likely to determine the saturation 
strength of magnetic field. Such a strong magnetic field is 



^ The viscous dissipation timescale is very long tv 

. q8 o23/4 rr.2 



( Shapiro|2o"o0 >. 
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Figure 3. Evolution of NS in the case when the Bondi-Hoyle regime begins at the surface of the RSG primary for the models RSG25He 
(left panel), RSG20 (middle panel), and RSG25 (right panel). The lines show the accreted mass normalized to Mq (solid line), the 
angular momentum of NS normalized to critical value jc = I^c, where flc = 760 Hz (dot-dashed line), and the orbital radius normalized 
to its initial value (dashed line). 




500 1000 1500 2000 2500 3000 
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Figure 4. The same as in Figure [3] for the case when the Bondi-Hoyle regime starts at a = 3 x 10 cm. 



also buoyant ( Spruitp999 ) and will emerge from under the 
NS surface into the surrounding accretion flow. 

Another possibility of generating super-strong magnetic 
in recycled NS was proposed by Spruit ( 1999 ) with appli- 
cation to X-ray binaries. In this model, the NS spins up 
to the critical rotation rate where it becomes unstable to 



exciting r-modes and gravitational radiation (Levin 1999 
Spruit|1999 ) . This may lead to effective braking of the outer 



layers of NS and developing of strong differential rotation 
Recent studies, however, suggest that non-linear coupling 
with other modes leads to the saturation of the r-mode at a 
much lower amplitude ( Arras et al.||2003 Brink et al.||2005 



Bondarescu et al. 2007), though the mass accretion rates 



assumed in these studies are much lower than those encoun- 
tered in common envelopes. 



3.2 Criteria for explosion 

Depending on the mass accretion rate, this may just modify 
the properties of the accretion flow near the star, the "accre- 
tor" regime, or drive an outflow (e.g. Hlarionov & Sunyaev' 



1975). In the latter case, two different regimes are normally 



discussed in the literature, the "propeller" regime and the 
"ejector" regime. In the ejector regime, the magnetosphere 
expands beyond the light cylinder and develops a pulsar 
wind. In the propeller regime, the magnetosphere remains 
confined within the light cylinder. If the magnetic axis is 
not aligned with the rotational axis it may act in a simi- 
lar fashion to an aircraft propeller, driving shock waves into 



the surrounding gas, heating it and spinning it up. If the 
energy supplied by the propeller is efficiently radiated away, 
this may lead to a quasi-static configuration (e.g. Mineshige 



et al. 1991). If not, the most likely outcome is an outflow. 



followed by expansion of the magnetosphere, and transition 
to the ejector regime. 

This classification is applied both in the case of quasi- 
spherical accretion and the case of disk accretion. In the disk 
case, the external gas above and below the disk is usually 
not considered as important dynamical element. However in 
our case this is not true because the accretion disk is rather 
compact and feeds from a more or less spherical component 
of the accretion flow. If the expanding magnetosphere re- 
moves this component, the disk will no longer be supplied 
with mass and quickly disappear, no matter what its regime 
is. For this reason, we will not consider the complicated in- 
teraction between the magnetosphere and the dislj^ but will 
focus on the interaction with the low angular momentum 
flow of the polar regions which accretes directly onto the NS 
avoiding the disk. In order to simplify the problem, we will 
treat this flow as spherically symmetric, keeping in mind 
that its mass accretion rate is only a fraction of the Bondi- 
Hoyle rate. 

The extent of the magnetosphere is often estimated by 
the balance between the magnetic pressure and the ram 



An interested reader may find out more about the disk- 
magnetosphere interaction from [Ustyugova et al.| ( [2b06| and ref- 
erences therein. 
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pressure of accreting flow. This is fine when the accretion 
shock is close to the NS. If not, the thermal pressure at the 
magnetospheric radius can be significantly higher compared 
both to the actual ram pressure at the shock and to the ram 
pressure that would be found at the magnetospheric radius 
if the free-fall flow could continue down to this radius. In 
this case, the initial radius of emerging magnetosphere is 
rather determined by the balance of the magnetic pressure 
and the thermal pressure of the quasi-hydrostatic "settling" 
flow downstream of the accretion shock ( |Mineshige et al | 
T99T|. 



Assuming that the accretion flow in the polar region is 
approximately spherical, we can describe it using the ana- 
lytical model developed by Houck Chevalier] 1991 ). Re- 



peating their calculations with the more accurate cooling 
function, adopted later in Chevalier (1996), 



10^ 



kT 
MeV 



erg 



(22) 



( Dicus|1972 Brown Weingartner|1994 ) and using the adi- 
abatic index 7 = 4/3, we find the accretion shock radius 



8.2x10" /ii?NsTMo"^/^'Mo"---cm, 



-1/27. 



-10/27 



where 



2R, 

Rns 



-1/2 



-64/27 



(23) 



(24) 



and Rg = GM/c? is the gravitational radius of the NS (In 
the calculations we assume the Schwarzschild space-time.). 
For reasonable masses and radii of NSs, 0.27 < fi < 0.62. 
In particular /i — 0.54 for M — Mq and i?Ns — 10 km. 
Downstream of the shock we have an adiabatic subsonic flow 
with the mass density 

-3 



where 



R 

Rsh 



M 



47r (GMy/^R: 



3/2 



is the mass density at the shock, the pressure 

-4 



P 



Pram(i?sl 



7 

where 

Pram(i?sh) = 



R 

Rsh 



47tr: 



,5/2 



(25) 



(26) 



(27) 



(28) 



is the ram pressure of the free falling flow at the shock, the 
temperature 



12 

11^^ 



1/4 



(29) 



where a is the radiation constant, and the radial component 
of 4- velocity 



M 



(30) 



Assuming that the magnetosphere expands until its 



magnetic pressure equals to the gas pressur^ we find the 
radius of dipolar magnetosphere with the surface strength 
B^s to be 

i?„ = 5 X 10Vr'^'i?N's(6 Mo-'/''Mo-'/''Bns,i5 cm. (31) 

For this to be above the NS radius the mass accretion rate 
should not exceed 

(32) 

yr 

The traditional criterion for the propeller regime is that 
the magnetospheric radius is below the light cylinder radius 



M„ - 1.3x10^ f-^'^'Bli%Rt,s,eMo 



i?Lc cP/27v ^ 4.8 X IO^P.3 cm 

and above the co-rotation radius 

1/3 



(33) 



Rc: 



fOMy 
V J 



1.5 X lO^P /^i?Ns,6 cm , (34) 



at which the magnetosphere rotates with the local Keplerian 
velocity. 

The ejector regime is realised when R^ > Rc. In terms 
of the mass accretion rate the propeller regime criterion can 
be written as 

and the ejector regime criterion as 



Mo 
yr 



M < Mej, 

where 
and 

M, ^ iT/r^^/^Mo-i?-^p^/^,,p-/^ ^. 

yr 

The latter is much stronger, and even if only a fraction of 
the total mass accretion rate is shared by the quasi-spherical 
component, it could be satisfied only during the early phase 
of the inspiral (see Figure [T]). More likely, the emerged mag- 
netosphere will operate in the propeller regime. 



(35) 



(36) 



The sound speed in the [Houck Chevalier (1991) so- 
lution depends only of the NS mass and the distance from 
the NS, 

c.^6.2xlO^Mo^/^i?-'/'— . (37) 
s 

The linear rotational velocity of the magnetosphere, 

^;rot ^4.8xl0^p-'i?6^, (38) 

is generally higher and thus we have the so-called super- 
sonic propeller. The power generated by such a propeller is 
roughly 



47r 2 

3 Cs 

s 



(39) 



^ As the result of its oscillations the NS can become quite hot 
and fill the magnetosphere with substantial amount of plasma. If 
this case the centrifugal force will need to be included in the force 
balance, yielding even higher value for R^. 
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( Mineshige et al.|1991| , exceeding by several orders of mag- 
nitude the accretion power 



GMM 



8xlO*''MoR^leMo—. 



(40) 



This suggests that the magnetosphere easily blows away the 
accreting matter and expands beyond the light cylinder, thus 
creating the conditions for proper pulsar wind from the NS. 

One may argue that the accretion flow could re- adjust 
to the new conditions, with the accretion shock moving fur- 
ther out, the temperature at the magnetospheric radius ris- 
ing, and the enhanced neutrino cooling compensating for 
the propeller heating, Ln = Lpr (A similar problem has 
been analysed by Mineshige et aL|(|1991 ).). Given the strong 
dependence of the neutrino emissivity on temperature, the 
cooling is expected to be significant only in the close vicinity 
of the magnetosphere. In the non-magnetic case, most of the 
cooling occurs in the volume comparable to the NS volume 



(Houck & Chevalier 1991). This suggests that in our case 



the cooling volume will be approximately the same as the 
volume of the magnetosphere. Then the energy balance can 
be written as 

|^i?ie„(T) = V, (41) 
where e^iT) is given by Eq 22 Solving this equation for R^, 



simultaneously with the pressure balance 



1 „2 ^ 



(42) 



(43) 



where p{R) is given by Eq 27 we find that 

-ORR^/l^ p-l/6 .^l/18p2/9 

-ft-NS 

independently on the mass accretion rate. Since in this solu- 
tion < Rns: we conclude that the accretion flow cannot 
re-adjust itself and will be terminated by the emerged mag- 
netosphere. 

In the case of cold neutron stars, the pulsar wind power 

s 



6xlO«BSs.i5flNs.eP-l'^,(44) 



where /x = BqR^q is the star magnetic moment and 0^ is 
the angle between the magnetic moment and the rotational 
axis of the NS (Spitkovsky"2006). For a hot NS the mass 
loading of the wind can be substantial, leading to even higher 
luminosity ( [Bucciantini et al.|2006l|Metzger et al.|2007| ). The 
energy released in the magnetically driven explosion is the 
rotational energy of the NS 



2 



1.6xlO''Moi?L,6J'-|erg, 



(45) 



where / ~ {2/5)M^sRls 

is the NS moment of inertia. 
Inside the RSG core the accretion rate is expected to 
be much higher. For example, in the RSG20He model the 
core density p ^ 8.6 x lO^g/cm^ and the pressure p ^ 1.3 x 
lO^^gcm^/s^, yielding the Bondi accretion rate 



Mb 



4.3x10^ 



yr 



(46) 



For such a high accretion rate, the propeller regime requires 
the magnetic field to be well in excess of lO^^G, similar to 
what is required in the magnetar models of GRBs. Whether 



such a strong dipolar magnetic field can be generated is not 
clear at present, as the observations of known magnetars 
suggest lower values. However, this could be a result of the 



short decay time expected for such a strong field (Heyl & 
|Kulkarni|1998l ). 



4 SUPERNOVA PROPERTIES 
4.1 The plateau. 

Since the NS has to reach the maximum allowed rotation 
rate before it explodes, the model predicts the explosion 
energy around lO^^erg, putting it on the level of hyper novae. 
The corresponding ejecta speed will be very high 

., = 10^M-y^HE, (47) 

where Mej,i is the ejected mass. 

Since the progenitor is a RSG the supernova spectra 
will show hydrogen lines and would be classified as type-II. 
As the NS spirals through the weakly bound outer envelope 
of its RSG companion, it drives an outflow with the speed 
about the local escape speed and the mass loss rate about 
Mq/jt (Taam & Sandquist 2000). The total in-spiral time 
is about few initial orbital periods. For a RSG with mass 
20 M© and radius R = lO^^cm this yields ^;esc = 5xl0^cm/s 
and tins ^ few x lO^s. Thus, we expect a hydrogen rich cir- 
cumstellar shell of few solar masses to exist at a distance of 
^ 10^^ cm at the time of the supernova explosion. The pres- 
ence of such a shell will delay the transition to adiabatic 
expansion phase and increase the supernova brightness sim- 
ilar to how this occurs type-IIn supernovae (Chugai et al. 



2004| [Woosley et al. 120071 [Smith et al.|2007[ |2008). 

Following Kasen & Bildsten (2010), the supernova lu- 



minosity during the phase of adiabatic expansion can be 
estimated as 



-E'rot te 



(48) 



where te = i?o / '^^ej is the expansion time scale and Rq is the 
progenitor radius, and 

1/2 



td 



3 Mej/C 

47r cVe\ 



(49) 



is the radiative diffusion time scale, k is the opacity (In the 
estimates below we use k — 0.4cm^/g, the value for Thom- 
son scattering in fully ionised plasma.). It is assumed that 
the initial magnetar spindown time scale is smaller com- 
pared to te, which is certainly satisfied in our case, given the 
large required mag netic field, 10^^ - lO^^G. For the charac- 
teristic parameters of the model we have 



te 

and 



10^i?o,i4^;e7,l s, td ^ 8 X lO'^Mll'^y-y^ s. 



1.6xl0^''i?0 14Mei 



erg 



(50) 



(51) 



This high luminosity will be sustained for up to t ^ td , after 
which the trapped radiation, generated by the blast wave 
when it was crossing the star, escapes and the luminosity 
drops. It most normal type-II supernovae the luminosity af- 
ter this point is sustained via radioactive decay of isotopes, 
mainly ^^Ni, produced by the supernova shock in the high 
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density environment around the collapsed core. However, in 
our model the stellar structure is significantly different from 
normal pre-supernovae. 



4.2 ^^Ni production 



If the explosion occurs after the NS is settled into the center 
of RSG20He the core density is only pc ^ 8xl0^gcm~^ and 
it temperature Tc 1.9 x 10^ K. where Ra = GM/c^ is the 
Bondi accretion radius. 

As the magnetically driven explosion develops a shock 
wave propagates through the body of RSG. Where the post- 
shock temperature exceeds Tni — lO^K, ^^Ni is produced in 



abundance ( [Woosley et al. 2002). In order to capture the 
shock dynamics, we need to know the structure of the ac- 
cretion flow onto the NS, which has been established before 
the explosion. This flow is rather complicated because of 
its rotation and without computer simulations we can only 
make very rough estimates. 

its rotation and near the NS there is an accretion disk. 
First consider the polar region where the flow is more of less 
spherically symmetric. Equation ] 13| shows that the accretion 
shock is quite close to the NS, particularly when it is near 
the RSG core, and hence most of the flow can be described 
using the free fall model. The ram pressure of the accretion 
flow in this model is given by 



2 MV2GM ^_5/2 

4:7T 



(52) 



where Va = {2GM / RY^'^ is the free fall speed. Close to the 
NS the speed of explosion shock v^i, is much less compared 
to and can be found from the pressure balance equation 



2 2 E 



7+1 



4nR^' 



(53) 



where R is the shock radius and E — L^t is the energy 
injected by the magnet ar. This yields 



t . 



2MV2GM 



(54) 



and thus the shock speed increases with the distance. This 
approximation breaks down at the radius 



Rt 



2MGM 



(55) 



where v^h = Va. The gas temperature exceeds Tni inside the 
radius i?Ni determined by the equation 



2 2 CLrr4. 



(56) 



and the total mass of produced ^^Ni can be estimated as 

MNi MtNi, (57) 

where t^i is the time required for the shock to reach i?Ni 
according to Eq|54] From these equations we find that 

Rt = 10^M6iv;;,5oCm, (58) 
i?Ni = 7.8 X lO^Me cm, (59) 

(60) 



MNi = O.OO3M0 M^^^^LZ!, 



Thus only a relatively small amount of ^^Ni can be produced 



even for the highest mass accretion rates allowed in this 
model. 

In the equatorial plane, the accretion shock can be 
pushed all the way towards R^cc (see Section[2]) . Downstream 
of the shock the flow is highly subsonic with 



(61) 



The highest density is achieved when the accretion shock is 
infinitely weak and this law is applied up to the sonic point 
Race- In this case, we can put Ro = i^acc and po = p*, the 
mass density of the RSG. The equation for shock speed is 
now 

2 2 ^ E 
7+1^"" 47ri?3' 

which yields the shock radius 



(62) 



R = 



2 f{7 + l)L^ 



1/2 



1/2 



3 V 87rp.i?Lc 
i?Ni is now determined by the equation 

2 2 aT^i 

where the mass of produced ^^Ni is still given by Eq, 
Combining these equations we find 



(63) 



(64) 



57 



(65) 



The highest values for Mni are obtained inside the RSG core. 
For the RSG20He, where i?acc 5xl0^cm, M - lO^Mo/yr 
(see Figures [l] and [2]), and ^ 5 x lO^g/cm"^, we find 
Mni — O.O2M0. Although significantly higher, this is still 
small compared to the mass deduced from the observations 
of most supernovae. 

If the density behind the shock wave is lower than 10^ g 
cm~^, which corresponds to M < lO^Mo/yr, then the pho- 
ton dissociation of Ni^^ becomes important as well ( |Magkot- 
|sios et a L]|2010|. Thus, our estimation gives only an upper 
limit on the amount of produced Ni^^. 

Another potential source of ^^Ni is the wind from the 
accretion disk itself, like in the collapsar model of GRB. 
Without a robust disk model, it is difficult to figure out how 
effective this mechanism is. For example, in the numerical 
simulations of [MacFadyen k, Woosley| ( |1999| ) , the outcome 
was dependent of the rate of introduced "viscous" dissipa- 
tion. In any case, the much lower accretion rates encountered 
in our model suggest that the amount of ^^Ni will be lower 
too. Thus, we conclude that our model predicts significant 
deficit of ^^Ni. 



4.3 The long term effects of magnet ar wind 

It has been suggested recently that the long-term activity of 
magnetars born during the core-collapse can power the su- 
pernova light curves, mimicking the effect of radioactive de- 



cay (Woosley 



2010 



Kasen & Bildsten 2010 ). In these studies. 



it was assumed that the rotational energy of the magnet ar, 
released at the rate oc t~^, is converted into heat at the base 
of the supernova ejecta. One mechanism of such heating is 
purely hydrodynamical, via the forward shock wave involved 
in the collision of the fast magnetar wind with the slower 
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supernova ejecta. In our case, however, this shock is the su- 
pernova shock and it is not present since the break-out. The 
other mechanism, heating by the radiation produced at the 
termination shock of the magnet ar wind, remains in force 
and here we explore its efficiency. 

Combining Eqs |44| and |45| it is easy to obtain the well 
known law of the asymptotic spin-down of magnetic rotators 

4- ) ^40B-l,,t-'^^s-' (66) 



and 



2fiH J 



-2 erg 



(67) 



(In this Section we use i?Ns = 10^ cm and Mns = 1.5Mq.) 
The last equation shows that unless B <^ 10^^ G the su- 
pernova luminosity will drop dramatically after the plateau 
phase. Such evolution has been observed in some super- 



novae, for example SN1994W (Tsvetkov] [1995| [Sollerman 
et al.||1998| ). 

In order to see how exactly this energy can be used to 
heat the supernova ejecta we can appeal to the observed 
properties of pulsar wind nebulae. For example, the study 



of the Crab Nebula by Kennel Coroniti| (.1984) indicates 
that immediately behind the termination shock most of the 
energy is in the form of ultra-relativistic electrons and only 
a fraction of it is in the magnetic field and there is no good 
reason to expect this to be different for magnet ar winds. 
However, the magnetic field of a one year old magnetar wind 
nebula (MWN) is expected to be very strong, making the 
synchrotron cooling time very short. 

One can estimate the magnetic field strength via the 
pressure balance at the termination shock 

where is the radius of termination shock. Assuming that 
— Vejt, this yields 



B 



0.3 tj ^ B^q iqV^^^q G . 



(69) 



The lowest energy electrons accelerated at the termi- 
nation s hock of the Crab Nebula have the Lorentz factor 
7 ~ 10^ ( Kennel Coroniti]|l984 ). In the magnetic field of 
the MWN, they produce synchrotron photons with energy 

^ 100 BoTiB^e^'/'t,-' keV , (70) 

and their synchrotron cooling time 

tsyn ^ 10^ t7^NS,15'^^ej,9 S , (71) 

is shorter compared to the dynamical time scale for few years 
after the explosion. 

The observations of the Crab Nebula also show that 
the electrons are accelerated up to the radiation reaction 



limit (de Jager et al. 1996 Komissarov & Lyutikov 2010). 
Assuming that the same is true in our case we expect the 
synchrotron spectrum to continue up to E^^^^^ ^ 100 MeV 
and the highest energy of the electrons to be 



100 trB 



1/2 1/2 



TeV. 



(72) 



These electrons also cool via the Inverse Compton (IC) 
scattering on the soft photons produced by the supernova 
shell. The ratio of the IC to the synchrotron energy losses 



^ = Usoit/uB, where Usoit is the energy density of the soft 
radiation and Ub — B^/87r is the magnetic energy density. 
During the plateau phase 



(73) 



where Rph is the radius of photosphere. Combining this re- 
sult with Eq. [68]we find 

2 



-^^soft / Rv 



(74) 



As we have seen, Lsoft can be much higher than Lw but on 
the other hand R^ can be much lower than Rph- Thus, it is 
hard to tell whether the IC emission can compete with the 
synchrotron one without detailed numerical simulation^ 
After the plateau phase i^soft drops dramatically and we ex- 
pect the synchrotron losses to be dominant. 

For photons with energy E^ < ruec^ - 0.5 MeV the 
main source of opacity is the Compton scattering in Thom- 
son regime. The corresponding optical depth is 



Tt ^ crTMej/47rmp(t^;ej)^ 



(75) 



where ctt is the Thomson cross-section. Thus, the shell be- 
comes transparent in about one year after the explosion. The 
typical energy of photons from radioactive decays is around 
1 MeV and hence the ejecta opacity to these photons is sim- 
ilar. As the result, the synchrotron photons and the photons 
from radioactive decays are utilised with similar efficiency. 

When El, ^ meC^ the scattering is in the Klein-Nishina 
regime and the cross section is smaller, ctrn ~ SaTmeC^ /SE^ 



( Lang|1980 ). For 100 MeV synchrotron photons this yields 

(76) 



Tkn 



10 ^Mej,lt7 '^ej,9 



and these photons begin to escape already after ten days or 
so. However, provided the spectrum of the shock accelerated 
electrons is a power-law N(E) oc E~^ with p > 2, which 
is expected to be the case, the contribution of such highly 
energetic photons to the energy transport is lower. 

The IC photons with energy E^ > {2meC?Y / E^^^oft — 
1.2 TeV will also interact with the soft supernova photons 
via the two photon pair-production reaction (here we used 
Ej^^soit = 1 eV). The corresponding opacity can be estimated 
as 

r^-f — ^ . . ^""T TP — — 2Lsoft,4i^^e7,U7^ • (77) 
5 47r(^;ejt)c^soft 

Thus, we expect the ejecta to become transparent to the IC 
emission soon after the end of the plateau phase. 

The above calculations show that during the plateau 
phase the high energy emission from the termination shock 
of the magnetar wind is deposited in the ejected. However 
this additional heating has little effect on the supernova lu- 
minosity. Indeed, the increase of the luminosity due to this 
heating can be estimated as 



AL '. 



(78) 



For t > td the energy deposition rate does no longer scale as 
(c.f. |Woosley,,2010J as the ejecta becomes transparent 



" Should the IC losses dominate, El 
reduced by the factors ^■'^/^ and ^ respectively. 



and E^^^^ are to be 
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to the high energy emission. For example, if the synchrotron 
emission is the main channel of heating the rate scales as 
LwT oc when r ^ 1, and the rate due to IC emission is 
likely to decline even faster. 

Out of the well studied supernova known to the au- 
thors none fits the above description. The closest example is 
SN1994W dCortini et al.|[T994l [SoTlerman et al^[T998| . This 
is one of the brightest type-II supernovae, with the peak 
luminosity ^ 10^^ erg/s. After around 120 days its luminos- 
ity drops dramatically down to ^ 10^^ erg/s. The tail of its 
light curve is very steep show ing very small mass of ejected 

In fact the data can be ap- 



Ni (Sollerman et al, 



1998) 



proximated by L oc t~ ' . All these properties agree nicely 
with our expectations. On the negative side, the spectral 
data indicate the lower eject a speed and the presence of a 
very slow component with v ^ 1000 km/s. This supernova is 
classified as type-IIn and seems to be explained very well in 
the model involving a collision between the SN ejecta and 
a massive circumstellar shell ejected few years earlier in gi- 
ant stellar eruption ( Chugai et al.|2004 Woosley et al.|2007 



Smith et al.||2007[ |2008t . 

The power of escaped high energy emission can be esti- 
mated as LwC"^. As usual in such cases, it will peak 
when T ~ 1, which may occur within the first year after the 
explosion. For a source at the distance d = 10 di Mpc, the 
corresponding total fiux will be 



(79) 



where tp is the time when the peak is reached. For the syn- 
chrotron component the spectral energy distribution is ex- 
pected to peak around = 100 keV (see Eq 70). 



4.4 Remnant 

What is left behind after the explosion depends on where 
inside the RSG it occurs. As we have demonstrated, this 
strongly depends on the mean specific angular momentum 
of gas captured via the Bondi-Hoyle mechanism. If it is as 



high as proposed in the early investigations (Illarionov & 



|Sunyaev||1975| [Shapiro Lightman||l976| then the NS will 
accrete only at the Eddington rate until it settles into the 
core of its companion. What would happens after this is not 
clear. The current view is that the neutrino cooling will be 
sufficiently high for the accretion to proceed at the Bondi 
rate and that the NS collapses into a black hole ( |Zel'dovich| 



et al. 1972). However, this conclusion is based on spheri- 



cally symmetric models. The high rotation rate developed 
in the core of the primary during the in-spiral means that 
instead of directly accreting onto the NS the core will form 
a torus, whose neutrino cooling may not be very effective, 
in which case such a conffguration can be relatively long- 
lived. Moreover, the mass accretion will be accompanied by 
efficient recycling of the NS, so one can still expect the NS 
to turn into a millisecond magnetar before its mass becomes 
high enough for gravitational collapse. Even if the magnet ar- 
driven explosion fails in the dense environment of the core, 
and the NS eventually collapses into a black hole, this will 
be a rapidly rotating black hole accompanied by a massive 
accretion disk - the configuration characteristic of the col- 
lapsar model for GRB. However, in contrast to the GRB 
progenitors the envelope of RSG star is not sufficiently com- 



pact for the relativistic jet to penetrate it during the typical 
life-time of the GRB central engine. Instead, the jet energy 
will be deposited in the RSG envelope and most likely re- 
sult in a supernova explosion, leaving behind a black hole 
remnant. Interestingly, the mass of the presupernova could 
be below the limit for black hole formation in the normal 
course of its evolution as a solitary star. 

If the angular momentum is much lower, as suggested 
bylDavies & Pringle (|1980| then the explosion occurs when 
the NS is still inside the RSG envelope then the remnant 
is likely to be a close binary consisting of a WR star and 
a magnetai[^ In order to show this, we first note that the 
pressure inside the core of RSG at the He-burning phase is 
Pc ^ lO^^ergcm"^ and much higher at the presupernova 
phase ( Woosley et al.||2002 ). This has to be compared with 
the pressure inside the high entropy bubble, created by the 
magnetar 

where Rh = Vht is the bubble radius, Vh is its expansion 
speed, and t is the time since the explosion. We are inter- 
ested in its value at the time when the bubble radius is com- 
parable with the separation between the NS and the RSG 
core. Taking Rp — 10^°i?b,io cm and t = Rh/vh we find that 



Ph • 



10''L^,49^^b,9 



qRu 



erg 



(81) 



If the explosion occurs while the NS is still well inside the 
RSG envelope then Rh,io ^ 1 and Pb ^ Pc- Thus, the 
magnetar wind cannot destroy the core. Another point is 
whether the gravitational coupling between the NS and the 
core is strong enough for the binary to survive the high 
mass loss during the explosion. Such a destruction is unlikely 
because for a sufficiently small separation at the time of the 
explosion most of the mass loss will come from the envelope, 
which has no effect on the gravitational attraction between 
the NS and the core. 

The magnetar wind will interact with the WR wind and 
will terminate when the radius of magnetar magnetosphere 
drops below the light cylinder radius. The magnetospheric 
radius is determined by the pressure balance 

^ ^^^^ ou^^^ 

(82) 



87r 



Rm 

Rns 



which yields 



Rrr 



; GxlO"" B 



1/3 



47ra2 



.V-l/6 -1/6 1/3 
5^*^WR,-5^WR,8^11 cm. 



(83) 



Using Eq|66] to find the light cylinder radius, we estimate 
the termination time of the magnetar wind as 



6w — ^NS,15^*^-^WR,-5^WR,8"ll • 



(84) 



The Bondi-Hoyle accretion radius of the magnetar in the 
WR wind 

2GM 



Ra — 



4xlO^°^;wR 8 cm 



(85) 



is above the magnetospheric radius, indicating that the mag- 
netar may begin to accrete when its co-rotation radius 



^ The 2D numerical simulations by [Fryer et al.| ( |l996| have in- 
dicated the possibility of a neutrino-driven explosion during the 
CE phase, also leaving behind a close binary system. 
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Rc^ 



(GM\ 



1/3 



dxlO"" Bli%tli'' cm 



(86) 



exceeds R^- Ignoring the decay of magnetar's magnetic field, 
we find that this will occur at 



t ^ 7x10^ B~l^^M. 



-1/2 -1/2 



WR,-5 ^WR 



,8«iiyr- 



(87) 



In fact, this time exceeds the expected decay time of the 
magnetar magnetic field 



lO'^.-'yr 



(Heyl & Kulkarni 1998). However, the relatively weak de- 
pendence of Rm on 5ns indicates that the effect of mag- 
netic field decay is minor and for the whole duration of the 
WR phase of the companion the magnetar will not be ac- 
creting. Instead, and will exhibit the behaviour typical for 
Soft Gamma Ray Repeaters or Anomalous X-Ray Pulsars 
(Thompson & Duncan 1996). 



Normally, a close high-mass binary system can produce 
at most two supernovae. However, our results suggest that 
in some cases this number can be increased to three. In- 
deed, in our model the first explosion creates the NS itself. 
When this NS spirals into its companion the second, now off- 
center explosion occurs. Finally, when the core of the WR 
star eventually collapsed, there is another explosion. More- 
over, the very close separation between the magnetar and 
WR star makes this binary system a promising GRB pro- 
genitor (e.g. Barkov Komissarov |[20T0 ) . Indeed, the WR 
star is very compact so the relativist ic GRB jets can break 
out in a reasonably short time and the spin-orbital inter- 
action ensures that the WR star is rapidly rotating. These 
are the two most important properties which a long GRB 
progenitor must have. 



5 CONCLUSIONS 

The main aim of this study was to investigate whether dur- 
ing the common envelope phase of a close binary, involving 
RSG and NS, the NS can spin up to a millisecond period 
and generate magnetar- strength magnetic field. If possible, 
this would result in magnetically driven stellar explosion, re- 
leasing up to 10^^ erg of magnetar's rotational energy inside 
the RSG. 

It turns out that the outcome is very sensitive to the 
specific angular momentum of the gas gravitationally cap- 
tured by the NS during the in-spiral. Should it be as high 
as suggested in the early papers by |Illarionov Sunyaev| 
(1975) and Shapiro & Lightman ( 1976| ), the accretion rate 
is low and neither mass nor spin of the NS will increase sig- 
nificantly until the NS settles into the center of RSG. This 
is because the accretion disk forms too far from the NS, its 
temperature remains rather low, and it is unable to cool 
effectively via neutrino emission, which is very sensitive to 
temperature. As the result, the accretion shock is pushed 
too far, beyond the sonic point of the Bondi-type flow, thus 
preventing the accretion from reaching the Bondi-Hoyle rate 
( Chevalier|1996t . The CE can either be ejected, leaving be- 
hind a close NS-WR binary or survive, leading to a merger 
of the NS with the RSG core ( |Taam &: San dquist 2000). It 
is not clear what exactly would occur following the merger, 
as we cannot treat this case using our approach. The com- 
mon belief is that the NS will begin to accrete at the Bondi 



rate and collapse into a black hole. However, this conclusion 
is based on models with spherical symmetry, whereas the 
strong rotation developed by the system during the merger 
makes the assumption of spherical symmetry unsuitable. It 
may still be possible that the NS spins up to a millisecond 
period and becomes a magnetar before its mass reaches the 
limit of gravitational collapse. 

If, on the contrary, the angular momentum is much 
lower, more in line with the analysis of Da vies &: Pringle] 
(1980]), the accretion can begin to proceed at the Bondi- 
Hoyle rate when the NS is still inside the common envelope. 
Further investigation of the Bondi-Hoyle accretion via 3D 
numerical simulations are required to clarify this issue. 

Once the NS becomes a millisecond magnetar, its 
emerging magnetosphere interacts with the accretion flow. If 
this occurs while the NS is still inside the common envelope 
and the magnetic field is as strong as ^ 10 ^^G, the inter- 
action is most likely to begin in the propeller regime and 
then quickly proceed to the ejector regime, given the small 
radius of the light cylinder. Strong magnetar wind blows 
away the common envelope and drives an explosion whose 
energy is likely to exceed the standard 10^^ erg of normal su- 
pernovae. However, the core of the companion survives and 
the explosion leaves behind a very close binary, consisting of 
a magnetar and a WR star. In spite of the very close separa- 
tion and the very strong wind from the WR star, the NS is 
shielded from the wind by its magnetosphere and prevented 
from accreting. Later on, when the WR explodes this will 
be a third explosion produced in the system. Moreover, the 
rapid rotation and the compactness of this WR star show 
that this explosion can be accompanied by a gamma-ray 
burst (Barkov & Komissarov 2010). 

If the magnetar forms only after the merger with the 
core then a higher magnetic field, ~ lO^^G, is required to 
explode the star. In this case, the remnant is a solitary mag- 
netar. If, however, the magnet ar-driven explosion fails and 
the NS collapses into a black hole, this will be a rapidly ro- 
tating black with a massive accretion disk. Although this is 
exactly the configuration proposed in the collapsar model 
of gamma-ray bursts, the relativist ic jets will not be able 
to escape from the extended envelope of RSG and produce 
such a burst. Instead, they will deposit their energy inside 
the envelope and drive a type-II supernova. 

Given the very high rotational energy of a millisecond 
NS the supernova is expected to be very bright at the plateau 
phase and show very broad spectral lines. However, only a 
small amount of ^^Ni is expected to be produced in the ex- 
plosion due to the relatively low densities, in comparison 
to those reached during the normal core-collapse explosions. 
By the end of the plateau phase the power of the magnetar 
wind is also significantly reduced and the supernova bright- 
ness is expected to drop sharply. The combination of the 
rapidly declining power of magnetar wind and the increas- 
ing transparency of the supernova ejecta to the high energy 
emission from the wind termination shock result in steeper 
than normal light-curve tails. 

The unique property of supernovae produced in this way 
is the high energy synchrotron and inverse Compton emis- 
sion from the magnetar wind nebula. Unfortunately, the su- 
pernova ejecta does not become transparent to this emission 
until a hundred days after the explosion, by which time the 
NS rotation is already very slow and its wind is no longer 
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that powerful. The flux of gamma-ray emission is expected 
to be rather low and difficult to observable, unless the ex- 
plosion occurs in the Local Group of galaxies. 
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